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.  Conditions  at  gun  muzzle 
conditions  along  trajectory 


I .  INTRODUCTION 


Long  rod  kinetic  energy  penetrators  often  find  tactical  application 
in  the  particular  aerodynamic  environment  of  normal  temperature  and 
pressure  and  flat  fire  trajectory.  This  follows  from  the  requirement  of 
very  high  velocity,  which  implies  a  short  time  of  flight  and  small  gravity 
drop,  and  from  the  fact  that  the  many  tank  battlefields  of  the  world 
are  at  nominal  sea- level  altitudes.  This  circumstance  permits  some 
simplification  to  the  calculation  of  the  trajectory  of  proposed  project¬ 
iles  still  in  the  preliminary  stages  of  design  and  analysis.  For 
preliminary  design  purposes,  the  Mach  number  excursion  can  be  considered 
very  short  and  the  drag  coefficient  over  the  range  of  flight  can  be 
considered  linear.  The  retardation  then  appears  in  closed  form  solution 
whereby,  given  the  drag  coefficient  characteristic,  the  velocity  becomes 
an  explicit  function  of  range. 

II.  PROCEDURE 

In  the  free  flight  regime,  the  force  balance  along  the  axis  of  the 
zero  yaw  trajectory  gives* 

F  =  m  v 

■-JPV2ACD  (1) 

where 

F  =  axial  force  opposing  the  acceleration  of  the  projectile, 
v  =  the  acceleration  of  the  projectile, 
p  =  the  density  of  the  flight  medium, 
v  =  velocity  of  the  projectile, 

A  =  reference  area,  and 
CD=  drag  coefficient. 

The  differential  expression  becomes 


d v 
v 


P  A  Cp  ds 
__ 


(2) 


where 


dv 

"  dt  ' 


ds  , 
v  "  dt  ’  and 

s  =  distance  along  trajectory. 


With 


CD=  k  M  +  b  , 

where 

b  =  intercept  of  vs  M  characteristic, 
k  =  slope  of  Cp  vs  M  characteristic,  and 

dv  2  ^ 

“v  “  2~m  p  A  M  +  b)  ds  =  -j£  ,  since  M  is  proportional  to  v, 

or 


dM  1 

“M  =  ‘  2“m  P  A  (k  M  +  b)  ds 
Upon  integration  (Appendix  A)  , 
b 


M  = 


R  eQs  -  k 


(3) 


where 

M  =  Mach  number  along  trajectory, 
s  =  distance  along  trajectory  , 
kM  +b 

R  =  0 

K  M 
o 

Mach  number  at  muzzle,  and 


Q 


P  A  b 
2  m 
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With  "s"  the  argument,  M  can  be  determined  directly.  This  equation  is 
presented  in  HP-97  program  formulation  in  Appendix  B. 

An  example  will  illustrate  the  application.  Figure  1-a  is  an  outline 
of  a  typical  long  rod  penetrator  (flechette)  for  which  range  data  and 
measured  retardation  to  600  meters  are  available2.  Figure  1-b  shows  an 
idealized  model  of  this  projectile  with  a  simple  fin.  The  corresponding 
experimental  and  calculated^  drag  coefficients  are  presented  in  Figure  2 
where  the  range  values  represent  turbulent  flow  conditions  for  small  yaw 
flight.  The  viscous  contribution  to  the  calculated  curve  assumes  turbulent 
flow  friction  factors  but  is  posed  for  zero  yaw  and  presumes  a  linear 
drag  characteristic. 


III.  RESULTS  AND  CONCLUSIONS 

In  Reference  2,  the  experimental  range  velocities  are  presented  to 
600  meters  (335,000  cal)  with  extrapolation  to  1000  meters  (536,102  cal). 
These  data  are  transposed  and  values  based  on  the  linearized  character¬ 
istic  and  Eq.  (3)  via  Appendix  A  are  given  on  Figure  3a  and  3b  for 

=  4.24  and  Mq  =  4.06  respectively.  The  maximum  deviation  is  within  3% 
over  the  full  range. 

Although  discrete  measurements  for  similar  projectiles  are  not  avail¬ 
able  for  comparison,  similar  agreement  may  be  expected. 


2 

Maynard  Piddington,  "The  Aerodynamic  Characteristics  of  a  SPIW  Projectile ", 
BRL  Memorandum  Report  1594 ,  September  1964.  (AD  #355679) 

William  F.  Donovan  and  Bertram  B.  Grollman ,  "Procedure  for  Estimating 
Zero  Yaw  Drag  Coefficients  for  Long  Rod  Projectiles  at  Mach  Numbers 
from  2  to  5 ",  ARBRL  MR  02819,  March  1978.  (AD  #A054326) 
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CALIBER  NOMENCLATURE  discussed  in  APPENDIX 


CALIBER  NOMENCLATURE  discussed  in  APPENDIX  C 


Figure  2.  Dreg  Coefficient  Characteristics 
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Figure  3-b.  Velocity  vs  Range  for  M  =  4.06 
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APPENDIX  A 


INTEGRATION  SUMMARY  OF  EQUATION  (3) 


From  Equation  (2) 


—  =  -  -r-  P  A  Cn  ds 
v  2  m  D 


=  -  p  A  (kM  ♦  b)  ds. 


dM  1  . 

M  (kM  +  b)  '  '  2  m  p  A  dSi 


and  #34  in  Pierce's  Table  of  Integrals*  gives 

-t  [*  -■  [*&]'' 

Jmi  L<  Jsi 


which  is  easily  verified  by  differei/tion. 


k  M  +  b 
o 


pA  Sj  b 


k  Mj  +  b 


for  s  =  0. 
o 


By  transposition  and  substitution 


*R.O.  Pierce,  A  Short  Table  of  Integrate ,  Ginn  and  Company,  Boston, 
1929. 

This  integration  form  was  independently  presented  by  Mr.  James  Bradley 
of  LFD. 


mm  mt  fxukdl 


p  A  b  s. 


k  M  +  b  „ - 

o  2  m 

- M - e 


R  e^sl 


k  Mj  +  b 

K 


ib 

ST 


,  and 


M, 


R  eQsl  -  k 


APPENDIX  B 


H-P  PROGRAM  FOR  RETARDATION 

This  program  requires  the  physicals  of  the  projectile  and  the 
linear  approximation  of  the  Cq  vs  M  curve.  It  includes  a  stepping 
feature  within  a  loop  which  automatically  decrements  the  range  and 
repeats  the  calculation.  The  input  is  dimensioned  in  the  units  in 
common  usage  locally. 

Input 

Rj  Extreme  range  in  meters 
R2  Initial  Mach  number 
R3  Shaft  diameter  in  inches 
R4  Projectile  weight  in  pounds 
Rs  Intercept  of  CD  vs  M  curve 
R^  Slope  of  CD  vs  M  curve 
Ry  Conversion  constant  =  .000671 
Rg  Range  decrement 


Output 

Range  in  meters 
Initial  Mach  number 
Mach  number  at  specified  range 
Retardation  in  ft/sec/km 
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APPENDIX  C 


CALIBER  NOMENCLATURE 

Caliber  nomenclature  is  widely  used  in  aerodynamic  expression  as  a 
dimensional  convenience  to  compare  performance  parameters  of  geometric¬ 
ally  similar  models.  It  is  usually  referred  to  a  linear  scale  represent¬ 
ing  the  arithmetic  ratio  of  a  linear  dimension  to  an  arbitrary  standard  - 
most  often  the  body  diameter  at  the  forward  bourrelet  -  but  has  been 
employed  to  identify  volumes*.  Only  a  simple  extension  of  the  reasoning 
is  required  then  to  simultaneously  de-dimensionalize  the  "mass"  factor 
in  a  given  expression  and  deduce  a  normalized  system  of  mechanical  units 
which  permits  a  rational  comparison  of  the  dynamic  properties  of  even 
geometrically  dissimilar  elements  of  machinery.  Usually  the  context  of 
discussion  identifies  the  quantities  as  "mass  cal",  "inertia  cal" 

"length  cal",  etc.,  although  a  complete  lexicon  of  explicit  and  descriptive 
terms  is  available  for  this  purpose. 

For  this  report,  the  following  correlation  is  employed: 


Length  (cal) 


linear  dimension 
diametral  dimension 


Weight  (cal3) 


gravity  weight _ _ 

gravity  weight  of  unit  volume  of  water 


-  S.G. 


Mass  (cal2  sec2) 


S.G. _ 

gravity  acceleration 


Thus,  with  force  equal  to  mass  times  acceleration: 


(cal3)  -  (cal2  sec2)  (ffi?) 


*MacAl lister,  et  al.,  "A  Compendium  of  Ballistic  Properties  of  Projectiles 
of  Possible  Interest  in  Small  Arms",  BRL  Report  No.  1S32,  February 
1971,  (AD  #882117). 
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2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


